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Abstract 



A self-consistent random phase approximation (RPA) is proposed as an ef- 
fective Hamiltonian method in Light-Front Field Theory (LFFT). We apply 
the general idea to the light-front massive Schwinger model to obtain a new 
bound state equation and solve it numerically. 
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I. INTRODUCTION 



One of the most interesting aspects of Light-Front Field Theory (LFFT) |TJ is that it 
shares many features with non-relativistic field theory, due to the simple vacuum. It makes 
the Hamiltonian method sensible and thus LFFT becomes an attractive framework for bound 
state problems. To my knowledge, however, only very limited techniques which are common 
in non-relativistic many-body problems have been used in the LF context, e.g., the Tamm- 
Dancoff approximation [0 and a boson expansion method pi. In order to fully exploit the 
advantages of LFFT, it is desirable to consider the possibility of the application of other 
many-body techniques to LFFT. 

In this paper, we consider a kind of random phase approximation (RPA) applied to bound 
state equations. In the usual RPA, only one-loop vacuum polarization diagrams (rings) are 
summed up Q. Here we replace the vacuum polarization with the "meson" intermediate 
states obtained by the bound state equation, as we did in a previous paper [Q. Thus our 
bound state equation is a self-consistent one. 

The motivation for this approximation are two-fold: (1) The most popular method of 
solving the bound state problems in LFFT is Hamiltonian diagonalization based on the 
Fock representation. In the crudest approximation, a meson is described as a qq valence 
state. (In general, such a description is very good in LFFT.) However it becomes rapidly 
more difficult to include higher Fock states. This difficulty leads to the idea of effective 
Hamiltonians |||7| which include the effects of higher Fock states. A problem in this approach 
is that the effective Hamiltonians are obtained only perturbatively. In particular, it has not 
been considered to include infinitely many Fock states. (2) In a previous paper || on the 
theta vacuum in the massive Schwinger model, we have shown that the vacuum polarization 
effects are very important for the low-energy physics. It is therefore natural to consider such 
effects in the calculations of meson states even in the absence of theta. Another important 
point in the previous work is that the vacuum polarization is approximately represented as 
intermediate "meson" states obtained by the (two-body) bound state equation. 

In order to make the idea as concrete as possible, we consider the massive Schwinger 
model as an example. In the next section, we review the perturbative Tamm-Dancoff trans- 
formation |J and apply it to the massive Schwinger model. In Sec. 3, we obtain the self- 
consistent RPA bound state equation. In Sec. 4, we numerically solve the equation. The 
final section is devoted to the conclusion and discussions. 

II. PERTURBATIVE TAMM-DANCOFF TRANSFORMATION 

In order to make the bound state calculation tractable, it is necessary to make the Hamil- 
tonian more diagonal in particle number space. In a previous paper ||, we eliminate the 
particle-number-changing interactions perturbatively by a similarity transformation. Al- 
though it was originally proposed to resolve the sector-dependent counterterm problem in 
LF Tamm-Dancoff approximation ||, it is also useful in the massive Schwinger model, where 
no ultraviolet divergence occurs. 

Consider the Hamiltonian H of the form, 

H = H + XW + gV (2.1) 
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where (H Q ) mn = Q m 5 mn . The interaction W does not change the particle number while 
V does. The parameters A and g are to organize the perturbation theory. We will put 
A = g = 1 at the end. 

We eliminate the particle-number-changing interaction perturbatively by the following 
similarity transformation, 

H eff = e iR He~ iR (2.2) 

so that H e ff is required to be diagonal in particle number. The operator R is expanded in 
a power series of g and A, 

R = gRi + g 2 R 2 + gXTt H . (2.3) 

To the second order in g, we obtain the following effective Hamiltonian, 

H eff = H + g 2 AH, (2.4) 

(AH) mn = W L i + ?T^7r) V "* V *n- ( 2 - 5 ) 

Although this effective Hamiltonian is rather simple, it however gets very complicated 
beyond the second-order in perturbation theory. Furthermore, it is perturbative by construc- 
tion and we do not have enough control over non-perturbative dynamics. For example, we do 
not know how to treat the non-perturbative divergences which occur in the diagonalization, 
except for the resummation by the use of the running coupling constant. 

In the case of the massive Schwinger model, however, there are no ultraviolet divergences 
and, as is seen in the next section, we can generalize these interactions (at least partially) 
to include some non-perturbative effects. 

For the massive Schwinger model, the effective interactions AH in the qq sector may be 
diagrammatically represented in Fig. [I]. 



III. A SELF-CONSISTENT RPA BOUND STATE EQUATION FOR THE 

MASSIVE SCHWINGER MODEL 

Among these effective interactions, we are particularly interested in the Coulomb inter- 
action with vacuum polarization being included. In Ref. ||, we showed that the vacuum 
polarization effects are very important in low-energy physics. It is therefore natural to 
include this effect (i.e., charge shielding) in the bound state calculations. Note that this 
approximation is essentially equivalent to the so-called random phase approximation (RPA) 
in many-body problems. Here, however, we do a bit further. As we did in Ref. 0, we re- 
place the vacuum polarization with an approximate complete set of meson states. A simple 
calculation gives the following replacement, 

-(2ir)e 2 Vkfkl5(k' + l' - k - I)-*— 

[k — i 

-> -(27i)e 2 Vk'rkl5(k' + I' - k - I) — — (3.1) 



e 2 



1 ' ' 2vr (k 




-i 
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where gv( n ) — Jo dxi/j n (x) and ip n (x) is the wave function of the n-th "meson" with mass 
fi n . (See Ref. for the convention.) k and I (k' and I') are longitudinal momenta of the 
incoming (outgoing) fermion and antifermion respectively. See Fig. |]. Note that we have 
summed the geometric series of the vacuum polarization ||. 

With this replacement, we obtain the following bound state equation (in units of e/ v^)? 

ii^ix) = -^—^(x) + I' dx^x) + f 1 dy ^f ~ ^ S({i> n }, {/.„}, m, x, y), (3.2) 



x(l — x) Jo ' Jo (x — y) 

where 

S({j>n},{Vn},m,x,y) = ll + -J2\9v(ri)\ 2 - ( ^_ y)2m2 + - | f x _ y) 2 m - ) \ (3-3) 

I n V^nH —y + (l-x)(l-y) ) ) 

is the screening factor. Eq. ( p.2|) is our extension of the 't Hooft-Bergknoff equation [§] , 
Some remarks are in order: 

1. The screening factor S depends not only on the momentum transfer (x — y) but on 
the individual momenta of the fermion and the antifermion. The x-dependent terms 
do not scale in a simple way. It makes the analysis of this equation very complicated. 

2. The equation ( |3.2| ) is a self-consistent one, so that we should solve it iteratively as 
we do in the next section. It is important to note that, although we started with 
the two-body bound state equation ( |3.4j ), our bound state equation ( |3.2| ) contains the 
effects of infinitely many Fock states because of this self-consistency. 



3. Although the equation ( |3.2| ) is a self-consistent one, it is not the usual Hartree-Fock 
approximation in which only the wave function of the specific state is included. Our 
equation requires the information on all the "meson" states. It is not the usual RPA ei- 
ther, in which only the perturbative vacuum polarization diagrams (rings) are summed 
and the self-consistency is not required. 

4. < S < 1. Of course it physically means that the electric charge is screened due to 
the vacuum polarization. The effective charge reduces in a momentum dependent way. 

5. The screening factor is included not only for the usual Coulomb term but also the 
so-called "self-induced inertia" term, 

i 1 \ 

dy- f ^ (3.5) 

o (x-yyj 

This is necessary in order to ensure the positivity of the eigenvalue /i 2 . By multiplying 
ijj'(x) and integrating over x from to 1, we obtain 
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/i 2 f l dxWx)\ 2 = m 2 t dx l f {x)l \ + 
Jo Jo x[l — x) 







dxijj(x] 



+ )- f dx f dy^p- — S({ip n },{[i n },m,x,y). (3.6) 
2 Jo Jo (x — yy 

Since the screening factor S is positive, the eigenvalue /z 2 must be positive. If we did 
not include the screening factor in the self-induced inertia term, we would not have 
had the positivity of /i 2 . 

6. In the zero- momentum-transfer limit (x — y) — > 0, the screening factor becomes mo- 
mentum independent, i.e., 



S({^ n l{ f , n lm,x,y)^\l + Y l ^ X ^\ ■ (3-7) 



|2 



-1 



The factor J2 n \9v{n)\ 2 / fJ^ is the one we encountered in Ref. B. 



7. The wave function tp(x) for the 't Hooft-Bergknoff equation ( ft.4|) has a definite meaning 
that it is the two-body amplitude with respect to the free fermion and antifermion. 
The wave function ip(x) for our equation ( |3.2| ) does not have such a clear meaning 
because it contains the contribution from the higher Fock states. Roughly speaking, 
however, it could be considered as the two-body amplitude with respect to the dressed 
(constituent) fermion and antifermion. 



IV. NUMERICAL SOLUTION 

In this section, we numerically solve our bound state equation ( |3.2| ). Unfortunately, 
however, we are unable to find a useful set of basis functions [|TTJ] which enable us to calculate 
the Hamiltonian matrix elements analytically. Instead, we use a naive discretization of 
the equation. We do not claim that the numerical results we present in this section is 
accurate. Rather we emphasize that our equation gives a reasonably small correction to the 
corresponding solution of the 't Hooft-Bergknoff equation and therefore our equation could 
be a starting point of further improvements. 

The method is extremely simple. We discretize the range of x, [0, 1], into N + 1 pieces, 
and impose the boundary condition on the wave function, ip(0) = if)(l) = 0. Thus our 
equation becomes an N x N matrix eigenvalue equation. We first solve it with 5 = 1, and 
then use the result to calculate the screening factor for the next calculation. In this way, we 
iteratively solve the equation until the solution becomes self-consistent. 

Most of the calculations were done with iV = 2000. We stop the iteration when the 
difference of the subsequent wavefunctions becomes less than e = 1.0 x 10~ 12 , 

f dx \i) new {x) - i) id{x)\ 2 < e. (4.1) 
Jo 

A typical number of iterations is 4. In general, the convergence is slower for smaller fermion 
masses. 
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First we compare the results with those without vacuum polarization. Fig. [3] shows how 
the meson mass calculated iteratively with the new bound state equation deviates from 
the one obtained with the 't Hooft-Bergknoff equation. We see that there is systematic 
improvement, which is almost independent of how fine the equations are discretized. For 
small fermion masses m ^ 0.2, because of the singular behavior of the wavefunction at 
the edges, even the result without vacuum polarization depends on N, and therefore the 
results there should not be trusted. See Fig. HI for the copmarison with the result obtained 



by the basis function method |10| which is accurate for small fermion masses. We however 
emphasize that the modification due to vacuum polarization does not ruin the spectrum but 
leads to a reasonable correction. This is not so trivial because the modification is essentially 
nonlinear and depends all the eigenvalues as well as the eigenstates. This, together with the 



valence dominance in LF massive Schwinger model [IT], gives a consistent picture that the 
modification with vacuum polarization correctly reflects the (small) effects of higher Fock 
states. 

Next we compare the wavef unctions. As seen in Fig. [5], the meson wavefunction for the 
new equation is less singular at the edges (x ~ and x ~ 1). 

As we emphasized earlier j8|, apparently we should not sum up the geometrical series 
of vacuum polarization. Here we did it as a regularization. In the calculations without 
summing up the geometrical series, we numerically found that a state whose wavefunction 
is odd under x <-> 1 — x becomes the lightest one when the fermion mass is small. This seems 
to be because 

1 ~~ o ^2 \9v(n)\ ( 2 ■ ( x _ y )2 m 2 H 2 : {x- y y7^~ ] (^-2) 

n \^n^ w A*n + (l-x)(l-y) I 

becomes very small when m — > 0. We are unable to give any positive interpretation of such 
a lightest state with wavefunction being odd under x <-> 1 — x, and thus think that it is an 
artefact. In fact, the odd lightest state occurs for m ^$ 0.1, where the results should not 
be trusted because, as we remarked earlier, the discretization cannot represent the correct 
endpoint behavior of the wave function, which is very important for the correct spectrum. 

We emphasized that the new bound state equation contains the information on all the 
"meson" states. It is however dominated by the lowest meson state and effectively reduces 
to a much simpler one. Namely, we numerically find that the sum in the screening factor 
3D is strongly dominated by the first term. This fact reduces the CPU time considerably. 
In Table |, we compare the reduced calculation with the full one. 



V. SUMMARY AND DISCUSSIONS 

In this paper, we developed a self-consistent RPA as a new effective Hamiltonian method 
in LFFT. The method was applied to the massive Schwinger model and a generalization of 
the 't Hooft-Bergknoff equation was obtained, which effectively includes the contributions 
from infinitely many Fock states. We numerically solved it and found that the modification, 
i.e., the inclusion of vacuum polarization, leads to reasonable small corrections. We hope that 
this new bound state equation becomes a starting point of more elaborate approximations. 

In the following, we would like to discuss several aspects about the new bound state 
equation. 
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It has been shown that the first-oder coefficient M\ of the expansion of the meson mass 
squared in terms of the fermion mass, 



^ 2 (m) = l + Mim + M 2 m 2 H , (5.1) 



cannot be improved even if the variational space is systematically extended |12| . Later, 
Dalley suggested that the endpoint behavior of the wave function may change if the effects 
of higher Fock states correctly. (It has been assumed that the endpoint behavior is solely 



determined by the two-body sector |IU,I2~|.) Because the coefficient is closely related to 
the endpoint behavior, it is hoped that the so-called "2% discrepancy" will disappear when 
the correct endpoint behavior is taken into account. Since our new bound state equation 
contains the effects of higher Fock states, it is natural to ask if the Dalley's suggestion is 
correct quantitatively. For small fermion masses, it looks legitimate to approximate the 
screening factor in the following way, 

W»},K},m,x,y)-5=|l + J^fflJ , (5.2) 

where fix = // is the meson mass and gy(l) = f dxtjj(x), and the bound state equation 
becomes 



= ~7^--M^ + T **K*) + ^ f dy M? ^ y)) . (5.3) 
x[l — x) Jo Jo {x — y) z 



We estimate the coefficient M\ by using the ansatz 

il)(x) = x p (l-x) f) , (5.4) 

and S = ~(1 + 0(m)), and find that Mi = 2%/\/6 = 2.56510 which should be compared 
with the value for the 't Hooft-Bergknoff equation, 27r/y / 3 = 3.62760 and the 'correct' value, 
2e 7 = 3.56125. Although it is considerably too small, it is important that it shifts from the 
value of the 't Hooft-Bergknoff equation. It is conceivable that the approximation employed 
here is too naive and a more elaborate calculation would give a much better value. In fact, 
the above approximation (the limit ( |5.2|) ) is NOT legitimate, because for finite m, however 
tiny, the term with the fermion mass is not negligible for x < 1/m 2 . At present, we are 
unable to extract any useful information on the endpoint behavior because of this subtlety. 
The numerical results in the previous section may suggest that the effects of higher Fock 
states would solve the 2% discrepancy, though the present numerical method is too crude 
to correctly calculate the meson mass in the small fermion mass region. 

Although we included the effects of vacuum polarization as an effective interaction, from 
the similarity transformation point of view, there are other kind of effective interactions. To 
ensure the invariance of the eigenvalues ("meson" masses) under the addition of effective 
interactions, it is necessary for the transformation to be unitary. That is, we have to add 
other interactions in Fig. |T] too. In this paper, we assumed that the effects of the other 
interactions are small. But this is a cheat: the self-energy diagram with vacuum polarization 
is infrared divergent. Presumably, such infrared divergences cancel in the calculation of 



the meson spectrum. (Some authors also assumed the cancellation [14||.) We are however 
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unable to pin point how such cancellation occurs. Note that the cancellation of infrared 
divergences would be different from that in scattering processes. In QED, for example, 
because positronium is a neutral particle and has a finite size, it cannot couple to very 
soft photons, which play an important role in the cancellation of infrared divergences in 
scattering processes. 

Finally, we give some comments on the outlook. 

1. It is of course interesting to extend the present formulation to realistic four-dimensional 
models. We are now working on this direction. 

2. It is interesting to note again that the new equation is a self-consistent one from the 
spontaneous-symmetry-breaking point of view. We hope that a similar setting would 
leads to a useful framework for the vacuum problem in LFFT. 

3. The vacuum polarization inside the meson is also important in the presence of theta. 
It solves several puzzles concerning theta. The work is now in progress. 
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FIG. 1. Effective interactions in the massive Schwinger model in the qq sector, generated by 
a perturbative Tamm-Dancoff transformation to order g 2 . Horizontal dashed lines stand for the 
instantaneous Coulomb interactions. 



10 



k' 



r 



z 



1 

FIG. 2. Effective Coulomb interactions with vacuum polarization being included. The geomet- 
rical sum of the vacuum polarization as an approximate complete set of "meson" states is taken. 
This figure is for k' > k and I > There is a similar (North-East) diagram for k > k' and /' > I. 
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Deviation from the ones without vacuum polarization 
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FIG. 3. The deviation of the meson mass with the new bound state equation from that with- 
out vacuum polarization. Data are for N = 2000 and for N = 3000. The deviation is almost 
independent of N. 



Limitation of the naive discretization 
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FIG. 4. The limitation of the naive discretization. The results with the naive discretization are 
compared with those obtained by basis function calculation for the 't Hooft-Bergknoff equation, 
the latter being considered to be accurate for small fermion masses. The naive discretization does 
not produce accurate values for small fermion masses. 
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wave functions for m=0.2 and for m=1 
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FIG. 5. The comparison of the wavefunctions obtained for the ("old") 't Hooft-Bergknoff equa- 
tion and those for the "new" equation for m = 0.2 and m = 1.0. The wavefunctions for the new 
equation are less singular at the edges for small fermion masses. 
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TABLES 



TABLE I. Comparison of the reduced calculation with the full one. The mass of the lowest 
meson state is shown. In the reduced calculation, we take only the first term in the sum in the 



screening factor ( |3.3| ). The calculations are done with the naive discretization explained in Sec. |V 
with N = 2000. 



fermion mass 


reduced calculation 


full calculation 


0.1 


1.07407458 


1.07407437 


0.5 


1.84809454 


1.84775659 


1.0 


2.82537904 


1.84775659 


5.0 


1.06492610 


1.06485924 
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